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A frequency domain model of multimode circular microresonators for filter applications in 
integrated optics is investigated. Analytical basis modes of 2D bent waveguides or curved 
interfaces are combined with modes of straight channels in a spatial coupled mode the- 
ory framework. Free of fitting parameters, the model allows to predict quite efficiently the 
spectral response of the microresonators. It turns out to be sufficient to take only a few 
dominant cavity modes into account. Comparisons of these simulations with computa- 
tionally more expensive rigorous numerical calculations show a satisfactory agreement. 

Introduction 

Nowadays due to their superior selectivity, compactness, and possibility of dense inte- 
gration, microresonators (MRs) become attractive for application as wavelength add/drop 
filters Q- A typical microresonator setting, where a ring/disk shaped cavity is placed be- 
tween two straight waveguides, is shown in Fig.[T] In this paper, we outline a spatial cou- 
pled mode theory (CMT) based model of 2D circular multimode optical microresonators. 




Figure 1 : Functional decomposition of 2D micro- 
resonators: R is the radius of the cavity with a core 
refractive index n c and a width w c (= 0, for disks). 
The straight waveguides have a core refractive in- 
dex n s and a width w s ; g\ and g2 are the separation 
distances between the cavity and the straight wave- 
guides; rij, is the background refractive index; L\ 
and Li denote the lengths of the cavity segments 
which are not included in the coupler regions. The 
letters A,B,C,D and a,b,c,d denote the coupler 
port planes. 



Abstract resonator model 



For modeling purposes, the MR is decomposed into two bent-straight waveguide couplers, 
represented by the blocks (I) and (II), which are connected to each other by two segments 
of the cavity. The external connections are provided by straight waveguides. Here we 
consider only forward propagating modes of uniform polarization. We assume that all 
elements are linear and that the back reflections inside the couplers and the cavity are 
negligible. Outside the couplers, the interaction between the constituent waveguides is 
assumed to vanish. We consider a frequency domain description of the optical field. The 
vacuum wavelength X prescribes the real angular frequency G). Assume that N s modes of 
the straight waveguides and Nb bend modes of the cavity are taken into account. 



Let y p be the complex propagation constant of the p'th cavity mode. The variables 
A q ,B q ,C q ,D q and aP ,b p ,c p ,d p denote the directional amplitudes of these properly nor- 
malized modes in the coupler port planes, combined into amplitude vectors A, B, C, D and 
a,b,c,d. Let S 1 and S n be the scattering matrices for coupler I and II respectively, i.e. 



S 1 



The amplitudes of the connecting cavity segments are related to each other as 



: p = b p e (-h"L l 
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Given input powers ff = \A q \ 2 at A and P\ = \C q \ 2 at C , we are interested in the trans- 
mitted powers Pj = \B q \ 2 at B and the backward dropped powers P^ = \D q \ 2 at D. This 
means solving the linear system of equations ©, for B q and D q , for q = 1,2, . . . ,N S . 
When scanned over a wavelength range, resonances appear as maxima of the dropped 
power and minima of the transmitted power. 

To evaluate the microresonator model described above, one must know the cavity prop- 
agation constants y p and the scattering matrices S 1 , S n . Using an analytic model of bent 
waveguides as described in Ref. Q, we obtain the bend modes and their propagation con- 
stants. Having access to the bend modes, a model of the bent-straight waveguide couplers 
in terms of CMT leads to the scattering matrices. A detailed description of this procedure 
for one cavity mode and one straight waveguide mode is presented in Ref. [3 1. In the next 
sections, we extend it to the case of multimode MRs. 



Multimode bent-straight waveguide couplers 

Consider the bent-straight waveguide coupler as shown in Fig. EKl). Let {E^,H^,£^} 
and {Ef ,Hf,£j} represent the modal electric fields, magnetic fields, and the spatial dis- 
tributions of the relative permittivity of the bent waveguide and the straight waveguide 
respectively. Here the modal fields include the harmonic dependence on the propagation 
coordinate and are expressed in the Cartesian coordinates (x,z)- The field {E,H} inside 

Figure 2: CMT setting for bent-straight 
waveguide couplers. The coupler is defined in 
the region [xi,x r ] x [z,-,z ]. The external seg- 
ments of the bent waveguide and the straight 
waveguide constitute the port connections. It is 
assumed that outside the coupler there is negli- 
gible interaction between the waveguides. 

the coupler is given by a linear combination of the modal fields of the bent waveguide 
(Fig.El-(2)) and the modal fields of the straight waveguide (Fig. El-(3)): 

E(x,z) 
H(x,z) 

where C?, Cf are unknown amplitudes. As shown in Ref. 0, a procedure based on the 
Lorentz reciprocity theorem leads to the following coupled mode equations 

M(x,z)-4C(z) = F(x,z)-C(z), (4) 
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+ E c q o 



Ef(x,z) 
Hf(x,z) 



(3) 



with My = (E^ , W m ; Ei , H£) = j a z • (E^ x Hf + Ejj* x H m )dx, Fy = / (e- £ m )E< 7 
C, = C'„ for /, j = 1,2, ■ • • , (A^ + A/j), and m = b if 1 < i < Nb otherwise m = s, n = b if 
1 < J < otherwise n = s. Here a z is a unit vector in z- direction and £ is the relative 
permittivity of the complete coupler. The integrations extend over [x[,x r ] for each z level. 

Solving Eq. © by a Runge Kutta method of order 4, we get a relation C(z () ) = T- C(z;) 
between the amplitudes in the output and input coupler ports. For rather radiative bend 
modes, it takes a long z- distance to stabilize the elements of matrix T. This difficulty 
is overcome by taking the projections of the coupled fields onto the straight waveguide 
modes ||4|. Then the output amplitudes of the straight waveguide modes are given by 



B q 



(Ef,Hf;Ef,Hf) 



(E?,H*;E?,H*) 
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(5) 



By incorporating these projection corrections into T, we finally obtain the required scat- 
tering matrix S. 



Simulations and comparison 

We consider a cavity in the form of a disk, i.e. w c = 0. Since the modal loss of the bend 
modes increases with growing mode order (where the order of the mode is defined as in 
Ref. [2|), one can expect that only the lower order bend modes play a dominant role for 
the field evolution in the cavity. 
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Figure 3 : TE power spectrum for a 2D disk microresonator. Left : CMT with only a single (one of TEo,TE i 
or TE2) cavity mode. Middle : CMT results for two (dashed line), three (solid line) and four (dash-dotted 
line) cavity modes together. Right: Comparison of FDTD results and CMT results with 3 cavity modes. 
MR specifications: n c = n s = 1.5, = 1.0, w c = (disk), w s = 0.4 /jm, R = 5 /jm, g = 0.2 /jm. 



The present CMT setting allows to investigate the significance of individual cavity modes 
for the spectral response of the MRs. The left plot of Fig. |3] shows the dropped power 
and the transmitted power when only a single cavity mode (either TEo, TEi or TE2) is 
included in the CMT model. The extrema corresponding to the fundamental mode (TEo) 
are much more pronounced than those related to the first order mode (TE1). If only the 
TE2 cavity mode is taken into account, on the present scale hardly any variations of Pd and 
Pj appear. The middle plot of Fig.|3]shows the cumulative effect of the higher order cavity 
modes. The extrema corresponding to the fundamental mode remain almost unaffected, 
but the shape of the resonances related to the TEi mode changes. Note that the curves for 
three and four cavity modes almost coincide. The right plot of Fig. |3] compares the CMT 
results (3 cavity modes) with FDTD simulations [5|. The results agree surprisingly well. 



Fig. IU illustrates the field profiles for the full MR structure as predicted by the CMT. 
At the resonance corresponding to the fundamental mode, most of the input power is 
coupled to the fundamental cavity mode and appears at the drop port. For the resonance 
corresponding to the higher order modes, the circular nodal line in the field pattern in the 
cavity indicates that a significant part of the input power is coupled to the TEi mode. For 
an off resonance wavelength, most of the input power appears at the through port. 
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Figure 4: Absolute value of the y component of TE fields of the MR. The gray scales (black = zero) are 
comparable among the plots. From left to right: MR field at a resonance corresponding to the fundamental 
mode (A, = 1 .043 /jm), a resonance corresponding to the higher modes (k = 1 .04833 p.m) and off resonance 
(k = l.055fjm). CMT simulations ( 3 cavity modes) for a MR with the parameters as for Fig. [3] 

Conclusions 

The CMT based model of 2D microresonators yields quite accurate results even if few 
cavity modes are taken into account. These results agree well with the rigorous FDTD re- 
sults and can be obtained with much lower computational effort. The role of the individual 
cavity modes can be clearly identified in this CMT model. 
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